Diffraction in Time: An Exactly Solvable Model 
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In the recent years, mater-wave interferometry has attracted growing attention due to its unique 
suitability for high-precision measurements and study of fundamental aspects of quantum theory. 
Diffraction and interference of matter waves can be observed not only at a spatial aperture (such as 
a screen edge, slit, or grating), but also at a time-domain aperture (such as an absorbing barrier, or 
"shutter", that is being periodically switched on and off). The wave phenomenon of the latter type 
is commonly referred to as "diffraction in time". Here, we introduce a versatile, exactly solvable 
model of diffraction in time. It describes time-evolution of an arbitrary initial quantum state in the 
presence of a time-dependent absorbing barrier, governed by an arbitrary aperture function. Our 
results enable a quantitative description of diffraction and interference patterns in a large variety 
of setups, and may be used to devise new diffraction and interference experiments with atoms and 
molecules. 

PACS numbers: 03.75.-b, 03.65.Nk, 42.25.Fx 



The wave nature of matter is both captivating and per- 
plexing, and its exploration has been at the center of 
experimental and theoretical research since early days of 
quantum mechanics. To date, diffraction and interference 
experiments have been successfully performed with parti- 
cles ranging from elementary particles and simple atoms 
to complex molecular clusters [HE]- Most spectacularly, 
wave-like behavior has been convincingly demonstrated 
even for such large organic compounds as C 6 o[Ci 2 F 2 5]io 
and C168H94F152O8N4S4, each comprising 430 atoms [3J- 

In optics, diffraction is typically portrayed as deflection 
of light incident upon an obstacle with sharp boundaries, 
that can not be accounted for by reflection or refraction. 
Interestingly, quantum mechanics allows for an addi- 
tional, intrinsically time-dependent manifestation of the 
phenomenon: Owing to the dispersive nature of quantum 
matter waves, sudden changes in boundary conditions 
may cause the particle wave function to develop interfer- 
ence fringes akin to those in stationary (optical) diffrac- 
tion problems. This phenomenon, pioneered in 1952 by 
Moshinsky [4 and presently referred to as "diffraction in 
time" (DIT), is at the heart of a vibrant area of experi- 
mental and theoretical research concerned with quantum 
transients (see Ref. [5] for a review). 

In the original Moshinsky's setup [I] , a monoenergetic 
beam of non-relativistic quantum particles is incident 
upon a perfectly absorbing barrier ( "shutter" ) . The shut- 
ter is suddenly removed at time t = t\. In other words, 
the transparency of the barrier x(i), called the aperture 
function, jumps instantaneously from at t < t\ to 1 
at t > t\, i.e., x(t) — @(* — £i)> with O denoting the 
Heaviside step function. The removal of the shutter ren- 
ders a quantum wave function with a sharp, discontin- 
uous wave front. The latter disperses in the course of 
time, smoothing out the initial discontinuity and devel- 



oping a sequence of interference fringes. As Moshinsky 
has shown, these fringes bear close mathematical similar- 
ity to those in Fresnel diffraction of light at the edge of a 
straight, semi-infinite screen. Moshinsky's paradigmatic 
triggered considerable interest to DIT in experimental 
research with ultra-cold neutrons and atoms, and Bose- 
Einstein condensates [H El E] • 

On the theoretical side, a number of interesting exten- 
sions and variations of Moshinsky's shutter problem have 
been addressed. Moshinsky himself extended his original 
result to the case of a "time slit" , in which the shutter 
stays open only during a time interval t\ < t < ti, as de- 
scribed by the aperture function x(t) = Q(t— ti)0(t2 — t) 
,6]. Scheitler and Kleber found an exact analytical solu- 
tion of a related problem, in which the role of a smoothly 
opening shutter was played by a time-dependent 5- 
potential barrier V(x,t) ~ t~ 1 S(x), with x denoting the 
spatial coordinate [318]. Various physical problems, re- 
quiring generalization of the original Moshinsky's setup 
to describe DIT caused by an arbitrary aperture func- 
tion, x(i), have been addressed in Refs. [SHU]- The 
analytical methods adopted in all these studies rely on 
treating the shutter as an effective "source" boundary 
of a semi- infinite coordinate space (transmission region) . 
The main advantage of this approach is that the trans- 
mitted wave can be readily expressed in terms of %(£) 
and a time-dependent source boundary condition. A sig- 
nificant difficulty with this approach however is that, in 
general, there is no unique well-defined recipe for finding 
the source function that would accurately mimic a given 
incident wave packet (although, in some cases, efficient 
approximations and exact bounds are known |12j). 

In this paper, we introduce a self-consistent exactly 
solvable model of DIT, free of arbitrary parameters. The 
model aims to describe dynamics of an arbitrary quan- 
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FIG. 1: Illustration of time-evolution of a quantum state in 
the presence of a time-dependent absorbing shutter. 



turn state in the presence of an absorbing time-dependent 
shutter characterized by an arbitrary aperture function 
x(t), see Fig. [I] It serves as a versatile generalization 
of the original Moshinsky's problem and reduces to the 
latter in the particular case x(t) = Q{t — ii). 

The central quantity analyzed in this paper is a prop- 
agator K(x,x']t) that relates the particle wave function 
$?(x,t) at time t > to an initial quantum state ^o(x) 
at t = through [T5] 
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The propagator satisfies the time-dependent Schrodinger 
equation (TDSE) on both sides of the absorbing shutter, 
positioned at x = 0, i.e., 



id t + dl) K(x,x';t)=0 for x,x'^0. (2) 

2 m 



Here, m denotes the mass of the particle. The propagator 
is subject to the initial condition 



K(x,x';0 + ) = 6(x-x') 



(3) 



and is required to vanish as x — > ±oo at negative imagi- 
nary times, i.e., at t = — i\t\. 

Our treatment of absorption is based on a time do- 
main rendition of an approach originally introduced by 
Kottler. The approach applies to diffraction of stationary 
fields governed by the Poisson's equation at spatial aper- 
tures in otherwise perfectly absorbing ("black") screens 
[TBI IT7] . According to this approach, a wave originating 
on one side of a black screen is subject to a discontinuous 
jump at every point of the screen, except for points inside 
possible openings (holes, slits, etc.). More concretely, 
the difference of the wave amplitude on the source side 
of the screen and that on the opposite side equals the 
amplitude of the corresponding free-space wave, i.e., the 
amplitude the wave would have if no screen were present. 
A similar discontinuity condition is imposed on the nor- 
mal derivative of the field. Both the field and its normal 
derivative are postulated continuous across the openings. 
As Kottler has shown, an exact solution of the Poisson's 
equation, subject to the above well-defined, though un- 
usual, discontinuous boundary conditions, is identical to 
the wave field predicted by Kirchhoff 's diffraction theory. 



In our model, we consider a time-dependent generaliza- 
tion of the Kottler's discontinuity condition at the shut- 
ter. Thus, we require the propagator to satisfy 

K(x,x';t)\lZt =sgn(x')[l- X (t)]K (x-x',t)\ x=0 (4) 
and 

d x K(x,x';t)\ X x Z - 

- sgn(x') [1 - X (t)]d x K (x - x', t)\ x=Q . (5) 

Here, x(t) denotes the aperture function, ranging be- 
tween (shutter closed) and 1 (shutter open), 
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is the free-particle propagator, and sgn(a/) = x'/|x'| 
stands for the sign function. The latter identifies whether 
the source is located to the left or right of the shutter. 

Equations (p]-[5]) completely specify dynamics of a 
quantum particle in the presence of a time-dependent 
absorbing barrier. In fact, as will be shown below, this 
quantum-mechanical problem admits an exact analytic 
solution valid for an arbitrary piecewise differentiable 
function x(0- The solution is given by 

K(x, x ; t) — S(x, x')Kq(x — x , t) + K\(x, x ; t) (7) 

with 

K 1 (x,x';t)= f dTuK (x,t-T)x(r)K (-x',T) (8) 
Jo 

sgn(.T') / x : 



and 



u(x, x';t,r) = 



(9) 



Here, E(x,x') — [l + sgn(x)sgn(iz/)]/2, i.e., S equals 1 if x 
and x' lie on the same side of the barrier, and otherwise. 
For what follows below, it is useful to rewrite the prop- 
agator, given by Eqs. ([7[|9]), in an alternative form. Eval- 
uating the integral in Eq. QSp by parts, we obtain 

K{x,x;t) = E(x,x')[l-x(t)]K a (x-x',t) + K 2 (x,x';t) , 

(10) 
where 



K 2 {x,x']t) 



x(o) + x(t) 
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and 
$(x,x';t,T) = 
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(12) 



Equations (10-12), while being mathematically equiv- 
alent to Eqs. Wj- 9h , enable a straightforward verification 
of the fact that K satisfies Eqs. (pHSl). Indeed, since 
both K 2 and d x K 2 are continuous functions of x, the va- 
lidity of the discontinuity conditions, Eqs. Q and (p)]), 
follows directly from Eq. (10). The validity of the ini- 
tial condition, Eq. p]), can be also verified straightfor- 
wardly: K(x, x'; 0+) = [S + (1 - E) X (0)}K (x, X 1 ; 0+) = 
S(x— x'). The fact that K satisfies the TDSE is confirmed 
by a direct substitution of Eqs. (10]-12) into Eq. (pi). 
(Here, it is convenient to take into account the identity 
9 t $ + ^a^ + i A-$(d x $) 2 = 0.) Finally, we note that 
uniqueness of the solution can be established in a stan- 
dard way [15] . 

It is important to point out that, due to nonunitarity 
of quantum evolution in the presence of absorption, the 
propagator K , in general, does not fulfill the composition 
property, i.e., 



+ OG 



K(x,x';t)^ d$K(x,&t-T)K(£,x';T). (13) 



This can be readily seen by considering the simple case 
of a time-independent, completely absorbing barrier, for 
which \ — and K = EKq. 

Although not valid generally, the composition property 
holds in the following important special case. Consider 
a scenario, in which the absorbing barrier acts only up 
to some time t{. In other words, suppose the aperture 
function \ is such that x(r) = 1 for all r > if. It can 
then be shown that K satisfies 



/oo 
d^K (x-^t-T)K^,x',r) 
-OO 



(14) 



for < t{ < t < t. The physical picture offered by 
Eq. (14) is apparent: Once the absorbing shutter has 
been switched off, the wave function evolves in accor- 
dance with the free-particle propagator. 

Equations ffiWty offer the following physical interpreta- 
tion of the wave function evolution. The full propagator 
if is a sum of Si^o, describing propagation in the case of 
a completely absorbing barrier (x — 0); arL d K\, repre- 
senting contribution of a barrier of nonzero transparency. 
In the transmission region, i.e., for x and x' such that 
E(x,x') = 0, the expression for K\, Eq. ^, conforms to 
the Huygens-Fresnel principle [18] : The probability am- 
plitude at the point x and time t, produced by a point 
source at x' and time 0, can be viewed as that produced 
by a fictitious source, located at the origin (between x 
and x'); the strength of the fictitious source is deter- 
mined by the free-particle wave function at the origin, 
modulated by the aperture function. 



In certain cases, the integral in Eq. ( 11 ) can be evalu- 



ated explicitly. One important example is that of a "time 
grating" , characterized by a "staircase" aperture function 

«+i - r) + 



XatO(t — tjy)i where < Xj < 1> with < j < N, and 
< U < £2 < • • • < iiV < t. In this case, the integral in 
Eq. (n]) reduces to X^=i(Xn+i - Xn)erf[$(cc,a; / ; t;t n )], 
rendering a closed form expression for the propagator. 
In particular, %(£) = Q(t — £1) leads to a propagator 
coincident with that in the original Moshinsky's shutter 
problem [T5] , 

A diffracted wave function, ^>(x,t), resulting from an 
arbitrary initial state, ^>q(x), can now be obtained by 
numerically evaluating the integral in Eq. (nl). As an 
illustration, we consider diffraction of an initially lo- 
calized Gaussian wave packet. Adopting atomic units, 
m = h = 1, we consider a coherent state ip qp {x) — 
(7r)~ 1//4 exp[— l(a; — q) 2 + ip(x — g)], with the average 
position q and momentum p, and set the initial state to 
be ^o( x ) = ip-w,5( x )- Figure [2] shows the wave function 
^(x,t) at time t = 3 diffracted at two different time grat- 
ings x( t )j shown in the insets of Figs. [2k and [2b. The 
elementary, one-period cell of the first grating, Fig. [2^,, is 
given by Ox 0(r)0(At-r) + l x 0(r- Ai)0(2A£-r), and 
that of the second grating, Fig. [2b, by 0x6(t)8(A£-t)+ 
0.5x8(r-At)e(2At-r)-l-lxe(r-2A£)e(3A£-r),with 
At = 0.056. (The value x = 0.5 means that the shutter 
allows only a half of the probability amplitude through, 
while absorbing the other half.) The probability densi- 
ties |\&(x, t)\ 2 for the two gratings are shown in the main 
panels of Figs. [2^, and [2b, and the corresponding Husimi 

distributions, H(q,p) 
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with 



the asterisk denoting complex conjugation, are presented 
in Figs. [2b and (2bl, respectively. (Phase-space repre- 
sentations, akin to H(q,p), have been previously used 
for establishing classical-quantum analogies in the dy- 
namics of diffracted particles [15].) The central bright 
peaks in Figs. [2t> and [2]i are localized around the point 
(q,p) — (5,5), which is the average phase-space location 
of the particle in the absence of a shutter. Additionally, 
in the case of the first grating, three separated diffrac- 
tion peaks - one in the reflection region, q < 0, and two 
in the transmission region, q > 0, - are clearly visible 
in Fig. [2Jd. On the contrary, in the case of the second 
grating, Fig. [2bl, diffraction peaks overlap substantially, 
forming a distinct pattern of probability naughts, i.e., 
points in phase space "avoided" by the diffracted parti- 
cle. We point out that, in general, the diffraction pattern 
sensitively depends on the initial wave packet and aper- 
ture function, making diffraction on time gratings well 
suited for spectroscopic analysis of quantum states. 

So far, we have shown that Eqs. ([7]-[9l, or equivalently 



Eqs. (|10|- 12|), give an exact solution to a quantum dynam- 

A natural question 



X(r) = Xo6(h - r) + E^i 1 XnO(r - t n )Q(t 



ical problem, defined by Eqs. ([2j[5[) 
is however in order: Does this dynamical system accu- 
rately model physical reality? In what follows, we argue 
that the proposed system should indeed provide a good 
description of the motion of a quantum particle in the 
presence of a purely absorbing (but not reflecting, in the 




FIG. 2: Diffraction in time of an initially Gaussian wave packet. Figures (a) and (b) show, respectively, the probability density 
and Husimi representation of the diffracted wave packet for the aperture function given in the inset of (a). Similarly, figures 
(c) and (d) correspond to the aperture function given in the inset of (c). All quantities are given in atomic units, m = h — 1. 



classical sense) barrier. The key building block of our 
model is a time-dependent extension of the Kottler's dis- 
continuity conditions, Eqs. Q and ^. The latter, in 
the case of stationary waves governed by the Poisson's 
equation, are mathematically equivalent to Kirchhoff's 
description of diffraction. Predictions of the Kirchhoff's 
theory, in turn, are generally found in good agreement 
with experimental data in the transmission region, while 
somewhat lacking accuracy in the reflection region |20j . 
Therefore, at the very least, it is reasonable to expect 
the propagator K(x, x';t) to correctly describe quantum 
wave dynamics in the transmission region. Moreover, if 
the absorbing shutter is in effect only until some final 
time if, so that x(t) = 1 f° r a U t > £f> then K(x,x';t), 
being an analytic function of x at any time t > if, must 
correctly describe the dynamics of the diffracted particle 
in both the transmission and the reflection regions. 

Nevertheless, in the absence of a compelling local the- 
ory of absorption, the ultimate answer to the question of 
how accurately the proposed model describes physical re- 
ality can only be given by an experimental investigation. 
To this end, atom-optics systems appear to be partic- 
ularly suitable. Remarkably, recent progress in control 



and manipulation of ultra-cold atoms has made it pos- 
sible to perform diffraction and interference experiments 
with a single, isolated atom, corresponding to a quantum 
wave packet highly localized in space [21]. At the same 
time, strong ionizing radiation has been used to realize 
an optical diffraction grating, similar in effect to an ab- 
sorbing nanostructure (mechanical) grating }22j . Utiliz- 
ing these technologies, one might envision a single-atom 
DIT experiment, in which quasi-one-dimensional motion 
of an atom is "intercepted" by a time-dependent absorb- 
ing shutter produced by a transversely-oriented, pulsing 
sheet of ionizing radiation. 

In conclusion, we have addressed a self-consistent 
mathematical model of diffraction in time and found its 
exact analytical solution in form of a time-dependent 
propagator. The latter enables a quantitative descrip- 
tion of diffraction of an arbitrary initial quantum state 
at an absorbing shutter governed by an arbitrary aper- 
ture function. We believe that our model will prove use- 
ful in the areas of coherent quantum control, quantum 
metrology, in designing new diffraction and interference 
experiments with atoms and molecules, and, more gen- 
erally, in exploring foundations of quantum physics. 



The author thanks Ilya Arakelyan, Orestis Georgiou, 
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